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We study the effects of diffuse Galactic, far-infrared extragalactic source, and radio 
point source emission on the cosmic microwave background (CMB) anisotropy data 
anticipated from the MAP experiment. We focus on the correlation function and genus 
statistics measured from mock MAP foreground-contaminated CMB anisotropy maps 
generated in a spatially-flat cosmological constant dominated cosmological model. 

^j- ' Analyses of the simulated MAP data at 90 GHz (0?3 FWHM resolution smoothed) 

show that foreground effects on the correlation function are small compared with cosmic 

£--. variance. However, the Galactic emission, even just from the region with \b\ > 20°, 

significantly affects the topology of CMB anisotropy, causing a negative genus shift 
non-Gaussianity signal. Given the expected level of cosmic variance, this effect can be 
effectively reduced by subtracting existing Galactic foreground emission models from 
the observed data. IRAS and DIRBE far- infrared extragalactic sources have little effect 
on the CMB anisotropy. Radio point sources raise the amplitude of the correlation 
I unci ion considerably on scales below 0?5, and also result in an asymmetry of the genus 
curve (at 2 a) on 0?5 scales. While the removal of bright sources above 5 a detection 
£> ■ limit effectively reduces the excess correlation, the genus asymmetry remains significant. 

Accurate radio point source data is essential for an unambiguous detection of CMB 
*-j \ anisotropy non-Gaussianity on these scales. 

When all the foreground sources are subtracted from the observed CMB anisotropy 
map, non-Gaussianity of cosmological origin can be detected at the 2 a level if the 
measured genus shift parameter |Az/| > 0.02 (0.04) or if the measured genus asymmetry 
parameter \Ag\ > 0.03 (0.08) on a 0?3 (1?0) FWHM scale. 

Subject headings: cosmic microwave background — cosmology: theory 
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1. Introduction 

Following the COBE-DMR experiment detection of anisotropy in the cosmic microwave back- 
ground (CMB) on large angular scales (Bennett et al. 1996; Gorski et al. 1996, 1998), there have 
been many measurements of CMB anisotropy on angular scales down to ~ 10' (see, e.g., Lee et al. 
2001; Netterfield et al. 2001; Halverson et al. 2001 for recent measurements of CMB anisotropy). 
These observations have begun to test cosmological models and provide interesting constraints on 
cosmological parameters (see, e.g., Kamionkowski & Kosowsky 1999; Rocha 1999; Gawiser & Silk 
2000 for reviews, and, e.g., Ratra et al. 1999a, 1999b; Rocha et al. 1999; Knox & Page 2000; Wang 
& Mathews 2000; Douspis et al. 2001; Podariu et al. 2001; Xu, Tegmark, & de Oliveira-Costa 
2001; Netterfield et al. 2001; Pryke et al. 2001; Stompor et al. 2001 for CMB anisotropy tests and 
constraints). 

A satellite mission, the Microwave Anisotropy Probe (MAP)^, will probe the CMB anisotropy 
at frequencies of 22, 30, 40, 60, and 90 GHz with expected angular resolutions of 0?93, 0?68, 0?47, 
0?35, and 0?21 FWHM, respectively. At 90 GHz the sensitivity of MAP is expected to be ~ 35 ^K 
per 0?3 x 0?3 pixel. The high quality CMB anisotropy data anticipated from MAP will significantly 
improve constraints on many cosmological parameters. 

Non-CMB foreground contamination must be understood and accounted for to fully utilize 
the cosmological discriminative power of MAP CMB anisotropy data. Recent rapid progress in 
the CMB anisotropy field has stimulated much work on modeling foregrounds and on developing 
methods for removing them from acquired data (see, e.g., Bennett et al. 1992, 1994; Brandt et al. 
1994; Dodelson & Stebbins 1994; Tegmark & Efstathiou 1996; Bouchet & Gispert 1999; Tegmark 
et al. 2000). There are two major potential sources for foreground contamination: diffuse Galactic 
emission and unresolved point sources. 

Kogut et al. (1996a, 1996b, hereafter K96) detect cross-correlations between the COBE-DMR 
and DIRBE maps, which they ascribe to Galactic dust and free-free emission. Using the Saskatoon 
data (~ 30 to 40 GHz, ~ 1° resolution, Netterfield et al. 1997) and the 19 GHz survey data (~ 3° 
resolution, Boughn et al. 1992), de Oliveira-Costa et al. (1997, 1998) 4 obtain a result similar to 
that of Kogut et al. Hamilton & Ganga (2001) find similar correlations using the UCSB South 
Pole 1994 data (~ 30 to 40 GHz, ~ 1° resolution, Gundersen et al. 1995; Ganga et al. 1997). At 
higher frequencies near 100 GHz and on similar angular scales, MAX data (Lim et al. 1996; Ganga 
et al. 1998) shows a similar correlation. On smaller angular scales Leitch et al. (2000) reach a 
similar conclusion using the OVRO data (~ 30 GHz). However, no such correlation is found in the 
Python V data (~ 40 GHz, ~ 1° resolution, Coble et al. 1999). Mukherjee et al. (2001) and de 
Oliveira-Costa et al. (2000) discuss foreground contamination of the lower frequency, larger scale 
Tenerife data. At frequencies higher than 60 GHz there is no significant correlation between CMB 
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maps and various synchrotron template maps such as the 408 MHz (Haslam et al. 1981) and the 
1420 MHz (Reich &: Reich 1988) radio survey maps (Kogut et al. 1996a; K96; de Oliveira-Costa et 
al. 1997, 1998). 

CMB anisotropy experiments with subdegree resolution may be sensitive to radio point sources. 
Gawiser &: Smoot (1997) study the effect of far-infrared (FIR) point source contamination of CMB 
anisotropy data by using the IRAS 1.2 Jy flux- limited galaxy catalog at various angular resolutions 
and observation frequencies. Gawiser, Jaffe, & Silk (1998a) derive upper and lower limits on 
microwave anisotropy from point sources over the range of frequencies 10 to 1000 GHz from recent 
sub-arcminute resolution point source observations. Their upper limit is ((AT/T) 2 ) 1 ' 2 ~ 10~ 5 for a 
10' beam at 100 GHz. Toffolatti et al. (1998) use a luminosity evolution model to make predictions 
for the number density distribution of extragalactic sources and their contribution to temperature 
fluctuations in each channel of the Planck Surveyor mission. Sokasian, Gawiser, & Smoot (1998) 
estimate the effect of bright radio point source contamination on MAP and Planck observations. 
They find that removing bright sources above the source detection limit (5 a — 1 Jy, where a is the 
rms level of CMB anisotropy measured) reduces the contamination from ((AT/T) 2 } 1 ' 2 ~ 10~ 5 to 
10~ 6 , consistent with the conclusion of Toffolatti et al. (1998). Recently, the Wavelength-Oriented 
Microwave Background Analysis Team (WOMBAT, Gawiser et al. 1998b; Jaffe et al. 1999) have 
presented many different all-sky high resolution CMB anisotropy contaminant foreground template 
maps, including those for Galactic dust, synchrotron, FIR point source, radio point source, and 
Sunyaev-Zel'dovich cluster emission. 

In this paper we generate mock all-sky foreground-contaminated CMB anisotropy maps for the 
90 GHz channel of MAP in a spatially- flat cosmological constant (A) dominated cold dark matter 
(CDM) cosmogony. We use these maps to investigate how foreground contamination affects the 
CMB anisotropy angular correlation function C(0) and the topology of CMB anisotropy. In §2 we 
summarize various foreground contamination sources. In §3 we summarize how mock foreground- 
contaminated CMB anisotropy maps are generated. In §4 we study the effects of foreground 
contamination on the amplitude and shape of the correlation function, focussing particularly on 
the location of the acoustic valley in the function 6 1 ' 2 C{6). In §5 we use the genus statistic to 
study the effect of foreground contamination on detecting non-Gaussianity of cosmological origin 
in the MAP data. We conclude in §6. 



2. Sources of Foreground Contamination 

2.1. Galactic Dust and Free-Free Emissions 

In the microwave sky three Galactic emission components, those due to synchrotron, free-free, 
and dust, are important. Although these Galactic foregrounds can be modeled theoretically pixel 
by pixel by multifrequency observations (Brandt et al. 1994), this requires very low instrument 
noise and accurate emission models. According to Bennett et al. (1992), synchrotron and free- 



free emissions have well-determined spectral behaviors with I v oc Z y^ - 7 ~^ 1 - 1 anc [ j u (x ^ _0 - 15 ) 
respectively. It is very difficult to produce an accurate all-sky dust emission model because its 
spectral behavior depends on the shape, composition, size distribution, and dynamics of the dust 
grains, and because the dust emission is a sum over the emission from all dust grains along each 
line of sight. 

Kogut et al. (1996a) and K96 apply a constant temperature warm dust and free- free emission 
model with spectral intensity 

Iu = r(u/u fB u (T) + A ff (u/u )-°- 15 (1) 

to the DIRBE maps and the DMR two- and four-year data. Here B U (T) is the Planck black body 
function, r opacity, (3 emissivity power-law index, vq = 900 GHz, and Aff the normalized amplitude 
of free-free emission. For the four-year DMR data they obtain r = (1.2+q;|) x 10~ 5 , T = 20.0110 
K, and f3 = 1.5l ' 3 (68% confidence limits) at \b\ > 20°. In eq. (1) the second term represents free- 
free emission, where the constant Aff is determined from the rms level of DMR-DIRBE correlated 
free- free emission, ATjj = 7.1 ± 1.7 /iK at 53 GHz (see §3 of K96). Since synchrotron emission 
is dominant only at frequencies lower than about 60 GHz, we neglect synchrotron emission in our 
analyses here. Figure 1 shows the spectra of the dust plus free-free emission model of K96. The 
thick solid curve represents the best-fit model. The upper and lower thin solid curves are the upper 
(T = 22.5 K, = 1.2) and lower (T = 16.0 K, (3 = 2.6) limits of the K96 model, r = 1.2 x 10~ 5 in 
all three cases. The dust temperature is almost inversely proportional to the emissivity (Kogut et 
al. 1996a). 

However, dust emission models with a constant dust temperature everywhere on the sky are a 
poor description of reality (e.g., Finkbeiner, Davis, & Schlegel 1999, hereafter FDS99). It is known 
from FIRAS observations that there seem to exist both warm (16 K < T < 21 K) and cold (4 
K < T < 7 K) dust components that contribute to Galactic dust emission (Wright et al. 1991; 
Reach et al. 1995). FDS99 provide a two-component temperature- varying dust model determined 
by fitting 100 //m dust emission data (Schlegel, Finkbeiner, & Davis 1998, described below) to 
FIRAS dust spectra. Their model parameter values are a\ = 1.67, «2 = 2.70, /i = 0.0363, and 
Q1/Q2 = 13.0, where a« is the i th dust component emissivity index, fi is the fraction of power 
absorbed and re-emitted by dust component i (with Ylkfk = -0' anc ^ 9* ^ s ^ e ra tio of the FIR 
emission cross-section to the UV/optical absorption cross-section for the i th component (see FDS99 
§4 for details). We use both the K96 one-component and the FDS99 two-component dust models 
in our analyses here. 

We need a high resolution Galactic emission map to model foreground contamination in the 
90 GHz MAP CMB anisotropy data. Schlegel et al. (1998) present a well-calibrated and Fourier- 
destriped full-sky 100 (iva. Galactic emission map constructed from the C05-B-DIRBE and IRAS- 
ISSA data sets, with zodiacal light and point source contamination removed. This map has an 
angular resolution of Q'A FWHM. We smooth this map using a Gaussian filter to get a map with 



total smoothing scale of 0?21 FWHM, and estimate the brightness temperature 5 at 90 GHz by 
using the Galactic emission spectrum template shown in Figure 1. For the FDS99 two-component 
dust model we estimate the 90 GHz dust emission by extrapolating the 100 (im flux using the dust 
temperature and DIRBE color-correction factor for each line of sight in the maps (see eq. [6] in 
§4 of FDS99). The averaged FDS99 dust emission spectrum (the thick long-dashed curve in Fig. 
1) is a factor of 1.4 lower than that of the K96 dust model. Since the FDS99 model does not 
include free-free emission, we generate a Galactic emission contamination map by adding free-free 
emission, assumed to be proportional to dust emission and to have an amplitude equal to 18% of 
the dust emission amplitude at 90 GHz, which is the percentage of the excess dust-correlated DMR 
emission estimated by FDS99 (see Table 4 of FDS99). Thermodynamic temperature fluctuations 
with respect to the mean CMB temperature To = 2.728 K (Fixsen et al. 1996) can be obtained 
from brightness temperature fluctuations multipied by the Planck factor (e x — l) 2 /x e x , where 
x = hv/kT®. This results in a high resolution all-sky temperature fluctuation map for Galactic 
emission that contributes to the MAP 90 GHz band data. 

While this is the best available map, it is not ideal. First, the K96 model parameters are 
determined from large (DMR) angular scale data while we study small-scale microwave temperature 
fluctuations. Second, the ratio of free-free to dust emission may be position dependent, resulting in 
a position-dependent Aff (while we assume that Aff is a constant). Third, if the free- free emission 
power spectrum is flatter than the dust emission power spectrum P(£) oc £~ 3 in multipole space, 
the ratio of free-free to dust emission will increase with decreasing angular scale (we assume that 
the large scale dust emission correlation with free- free emission persists to smaller scales). 

In this paper we use the K96 model free-free emission term to account for about 20% of 
the dust-correlated emission (see FDS99 for a related discussion). Similar correlations have been 
detected for many CMB observations (see §1, and Kogut 1999 for a review). This excess microwave 
emission is clearly correlated with the dust but is not due to thermal (vibrational) dust emission. 
Draine & Lazarian (1998) propose that K96's dust-correlated emission is not free- free emission but 
emission from rapidly spinning dust grains. Nevertheless, even if the dust-correlated emission is 
not free-free emission, our phenomenological approach is fairly reasonable, especially at 90 GHz. 



2.2. Far Infrared Extragalactic Sources 

The IRAS 1.2 Jy flux-limited galaxy catalog (Strauss et al. 1990; Fisher et al. 1995) contains 
5,319 galaxies, almost all of which are inactive spiral galaxies though some are active star-forming 
galaxies and active galactic nuclei. The IR emission from inactive spiral galaxies is from reradiation 
of interstellar dust heated by the general interstellar radiation field (Helou 1986). Assuming that 
these IRAS galaxies are all inactive spirals with emission spectrum similar to that of our Galaxy, 



'Brightness temperature Ta is related to intensity I v through I v — 2kTA(v 2 / (?) where k is Boltzmann's constant. 



we transform their 100 /xm fluxes to 90 GHz by using the K96 Galactic emission spectrum. We 
also use the method of Gawiser &; Smoot (1997) and Gawiser et al. (1998a) to determine the 
rms thermodynamic temperature fluctuation caused by IRAS extragalactic sources, and obtain 
((AT/T) 2 ) 1 / 2 ~ 5 x 10~ 7 at 90 GHz and 0?21 FWHM angular resolution. This is consistent with 
the result of Gawiser h Smoot (1997). 

A second set of contaminating FIR extragalactic sources are the few hundred sources we 
have tentatively detected in the DIRBE maps. DIRBE observed the whole sky with a 0?7 x 0?7 
square beam and produced infrared maps at 10 near-to-far infrared wavelengths, with pixel size of 
0?32 x 0?32. We use the ZSMA (Zodiacal light-Subtracted Mission Average) maps at 12, 25, 60, 
100, 140, and 240 //m to search for dim signals in the maps. We generate a difference map at each 
band by subtracting the intensity map smoothed over 2?36 FWHM from the original unsmoothed 
map. Thus only signals on scales smaller than 2?36 remain in the difference maps. The i pixel 
noise ctj in the difference map is Oi = (Tw/s/Ni, where &n is the global noise level of the difference 
map and TV, is the number of observations at the i th pixel, gn can be determined iteratively from 
the less-contaminated background region (|6| > 60° and ecliptic latitude \(5\ > 20°) by requiring 
z = lijoi (where li is the intensity at i th pixel) to have a normal distribution with zero mean and 
unit variance at z < 0. The resulting global noise levels are 0.209, 0.295, 0.549, 1.14, and 34.8, and 
19.9 MJy/sr at 12, 25, 60, 100, 140, and 240 //m, respectively. The noise levels in the 140 and 240 
^m difference maps are very large so we exclude these maps from our analysis. 

In the region |6| > 20° of the 100 //m map we select pixels with i^ioo /an > 3ctj. A faint 
extragalactic source may appear in at most 4 pixels (0?64 x 0?64) in the DIRBE map, when it is 
located where the 4 pixels touch. By grouping contiguous pixels, we search for signals in groups 
of 4 or fewer pixels in the Quadrilateralized Spherical Cube projection (see COBE DIRBE Exp. 
Supp. 1998). This excludes most of the bright galaxies detected by Odenwald, Newmark, & Smoot 
(1998), since they are typically detected in more than 4 pixels. We find 421 and 401 such groups 
(with 4 or fewer pixels) in the Galactic north and south hemispheres, respectively. 436 of these 
groups have positive fluxes at all wavelengths considered. We finally catalogue 289 extragalactic 
source candidates satisfying the IRAS bright galaxy color criteria (Smith et al. 1987; Soifer et al. 
1989; Strauss et al. 1990). These color criteria are —0.9 < log(/6o ^m//ioo ^m) < 0.3, —1.5 < 

l0g(/25 ina/fw ,um) < 0.0, -1.2 < log(/i 2 M m//25 /urn) < 0.4, -2.1 < log(/i 2 ^m/ZeO fim) < 0.0, 

and /| > /12 ^m/25 /urn- Figure 2 shows the distribution of the final 289 DIRBE FIR sources. 
33 IRAS galaxies and only one source in the 5 GHz 1 Jy flux- limited extragalactic radio source 
catalog (518 objects, Kiihr et al. 1981) are within a radius of 0?32 of these DIRBE sources. Most 
of these sources are matched with extragalactic objects in the IRAS Faint Source Catalog within 
the same radius. However, the large DIRBE pixel size prevents us from uniquely identifying these 
sources as extragalactic because each DIRBE pixel contains many IRAS sources, some of which 
are not extragalactic. Some of the DIRBE sources could conceivably be objects in our Galaxy 
that coincidentally satisfy the IRAS galaxy color criteria. The emission from these DIRBE FIR 
sources (after exclusion of the 33 sources identified with IRAS galaxies) are transformed to 90 GHz 



thermodynamic temperature fluctuations using the method applied above to the IRAS galaxies. 



2.3. Radio Point Sources 

The 4.85 GHz Green Bank survey catalog (Becker, White, & Edwards 1991) contains 53,522 
radio point sources with fluxes larger than about 40 mJy and covers 0° < 5 < 75°. The 4.85 GHz 
Parkes-MIT-NRAO survey catalog (Griffith k, Wright 1993) has 50,814 point sources with a flux 
limit of about 35 mJy in the — 87?5 < 5 < 10° region. Besides these, there is a 1.4 GHz survey 
catalog that has 30,239 sources flux-limited at 100 mJy in the —5° < 5 < 82° region (White & 
Becker 1992). These catalogs contain extragalactic sources as well as Galactic radio sources. For 
example, the Green Bank survey catalog has about 350 galaxies, ~ 800 quasars, more than 500 BL 
Lacertae objects, and ~ 400 X-ray sources (Becker et al. 1991). 

To make a combined radio source catalog that almost covers the sky, we adopt the Green 
Bank catalog for 0° < 5 < 75°, the Parkes-MIT-NRAO catalog for -87? 5 < 5 < 0°, and the 1.4 
GHz radio source catalog for 75° < 5 < 82°. In the few GHz frequency band, the main emission 
mechanism for radio sources is synchrotron radiation with power-law spectral behavior I u ~ v a with 
— 1.0 < a < —0.5 (Platania et al. 1998). Since spectral indices of most sources in the Green Bank 
and Parkes-MIT-NRAO catalogs are close to —0.85, radio sources with unknown a are assumed 
to have a = —0.85. Some radio sources with compact active nuclei generate flat-spectrum radio 
emission (a > 0) due to self-absorption of low energy photons. The average spectral index of flat- 
spectrum compact sources is a ~ at v < 200 GHz, steepening to a = —0.7 at higher frequencies 
(Danese et al. 1987; Toffolatti et al. 1998). In our analyses here we set a = for all radio sources 
with spectral index a > 0. Emission from radio sources in the region 75° < 5 < 82° are transformed 
to 4.85 GHz fluxes using the spectral index determined above. The fluxes of the final 72,132 4.85 
GHz sources, covering — 87? 5 < 5 < 82°, with flux limit of 43 mJy, are converted to temperature 
fluctuations when "observed" with the 0?21 FWHM MAP beam at 90 GHz. 



3. Generating Mock MAP Maps 

Park et al. (1998, hereafter P98) generate mock CMB anisotropy maps for the 90 GHz (0?21 
FWHM beamwidth) channel of the MAP experiment in a spatially-flat A dominated CDM model. 
This ACDM cosmogony has non-relativistic matter density parameter Qq = 0.4, h = 0.6 (where the 
Hubble constant Hq = 100/t km s _1 Mpc^ 1 ), baryonic-matter density parameter Qb = 0.0125/i -2 , 
and is normalized to the four- year DMR data (Ratra et al. 1997). The rms level of the CMB 
anisotropy in this map is 104 ^K (at 0?21 FWHM), larger than that of the flat CDM case (92 ^K; 
Refregier, Spergel, & Herbig 2000). The maps include MAP instrumental noise, simulated to give 
35 fiK per 0?3 x 0?3 pixel. 

It is expected that MAP will produce CMB temperature fluctuation maps with zero mean even 



if foreground contamination is not completely removed. Thus in the map the observed temperature 
at each position on the sky, T(0, (ft), is the sum of the CMB contribution Tq + STqmb, a contribution 
from foreground contamination STq, and an instrumental noise contribution <5T/v, where To is CMB 
black body temperature (2.728 K) and 5Tcmb is the CMB anisotropy temperature. Therefore the 
observed temperature anisotropy is 

ST T{0,4>)-{T) 5T C mb(0, 4>) + ST c (0, <p) + 5T N (0, 0) - (ST C ) 

jri.0, <t>) Jf-) " T (2) 

where (• • •) represents an average over the whole sky excluding the Galactic plane region (typically 
the \b\ > 20° part of the sky in our analyses here), and {STcmb) = (STpj) = are assumed. 
The term (STq) that will appear in the denominator of the last expression is omitted, since it is 
negligible compared with To. 

We use eq. (2) to generate foreground-contaminated CMB temperature fluctuation maps at 
90 GHz with 0?21 FWHM resolution. These maps are further smoothed over 0?214 (i.e., a total 
smoothing scale of 0?3 FWHM) to reduce noise, and consist of 3600 x 1801 points in (£, b) space 
with 0?1 spacing. Figure 3 shows one of the foreground-contaminated CMB temperature fluctuation 
maps (using the K96 Galactic emission model). The Galactic emission model has an uncertainty. 
Its uncertainty (hereafter AG) is taken to be the difference between the maps generated from 
the K96 (constant dust temperature) and FDS99 (spatially-varying dust temperature) Galactic 
emission models. Table 1 lists the map rms values of the CMB temperature fluctuations, MAP 
instrumental noise, and various foreground emissions. Note that rejecting radio sources above 5 a 
reduces the rms level of temperature fluctuations. When analyzing the maps we mask the LMC, 
SMC, M31 and Orion- Taurus regions for Galactic emission-contaminated maps. For radio source 
emission-contaminated maps we mask regions that are not surveyed (6 > 82°, 5 < —87.5°) or only 
partially surveyed. 



4. Correlation Function 

4.1. Amplitude of Correlation Function 

The two-point angular correlation function (hereafter CF) for temperature anisotropy is C{0) = 
(5T/T(6i, (pi)5T/T(6j 1 4>j)), where (• • •) denotes an ensemble average and is the angular separation 
between (0i,(pi) and (0j,(j)j). It is related to the angular power spectrum Cg through 

«2L Of 4- 1 
C ^ = E -^C e Pt(<x*6)BiFf, (3) 

1=2 

where Pi is a Legendre polynomial and Bg and Fg are the coefficients of the beam function and 
smoothing filter expanded in a Legendre polynomial basis. 

We measure C{0) directly from the CMB maps by picking up a million pairs of points for each 
angular bin considered. Figure 4o shows the angular CFs of noise-free and noise-added ACDM 



CMB anisotropy maps. The theoretical CMB anisotropy CF of the ACDM cosmological model 
(dot-dashed line) is also shown for comparison. Note that the amplitude of the CF can exhibit 
large cosmic variance while the position of the acoustic valley does not (see §3 of P98 for the effect 
of cosmic variance on the position of the acoustic valley of the CF). The amplitude of the CF of 
the noise-added CMB map is slightly higher than that of the noise-free case on scales less than 0?3. 
The effect of noise on the CF amplitude exceeds that of cosmic variance on these scales. On scales 
larger than 0?3, however, instrumental noise hardly affects the CF and cosmic variance dominates 
the statistical uncertainty. The effect on the CF of Galactic plane cuts at \b\ = 20° and 40° is 
illustrated in Figure 46. 

The effects of foreground emission on the CF are shown in Figure 5, where the CFs are 
computed for the Galactic plane cut region with \b\ > 20°. Although foreground Galactic emission, 
modeled either by K96 or FDS99, slightly increases the correlation amplitude on all angular scales 
(Fig. 5a), these effects are smaller than the cosmic and sample variance uncertainties. Figure 5b 
demonstrates that the Galactic emission model uncertainty AG (i.e., Galactic emission residuals 
which remain after removing model Galactic emission from a foreground-contaminated map) hardly 
affects the CF. IRAS and DIRBE FIR extragalactic sources appear not to have any effect (Fig. 
5c). Radio point source emission raises the CF amplitude on scales smaller than 0?5. When bright 
radio sources above 5 a are removed this small-scale effect becomes much weaker (Fig. 5d). 



4.2. Position of the Acoustic Valley 

The Doppler peaks in the angular power spectrum, particularly the first Doppler peak, appear 
as an acoustic valley in 6 1 ' 2 C(0). The position, depth, and width of the acoustic valley depend on 
the cosmological model parameters (P98) . In addition to the effect of foreground contamination on 
the amplitude of the CF, it is necessary to study their effect on the position of the acoustic valley. 

To predict acoustic valley positions we compute C{6) (using eq. [3]) from 200 sets of aJT'B^ 
and (a™ Be + n)Fn up to £ = 2500. We assume the noise n is a Gaussian random variable with 
zero mean and variance appropriate for MAP, and here Bi and Fg represent the 0?21 FWHM beam 
and 0?214 FWHM smoothing filter, respectively. The uncertainty in the position of the acoustic 
valley was estimated from these 200 Monte Carlo realizations (P98). The ACDM model predicts 
the location of the acoustic valley to be at V = 1?12 ± 0?02 for a 0?21 beam. When we account 
for the MAP noise and use an additional 0?214 FWHM smoothing filter (total smoothing of 0?3 
FWHM), the position becomes V = 1?14 ± 0?02 while it is 0?65 ± 0?04 for an open CDM model 
with Q = 0.4, h = 0.65. 

We have also computed C(9) near the acoustic valley directly from the foreground-contaminated 
CMB anisotropy maps, and have measured the position of the acoustic valley. The results are shown 
in Figure 6 and Table 2. The location of the acoustic valley is quite robust. When foreground Galac- 
tic emission is reasonably modeled, the Galactic cut mainly affects its location (~ 0?01). The cosmic 
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variance is about 0?02. 



5. Topology 

5.1. The Genus 

Gaussianity and phase correlation information are important probes of the origin of structure 
formation. Topology measures for two-dimensional fields, first introduced by Coles & Barrow 
(1987), Melott et al. (1989) and Gott et al. (1990), have been applied to observational data such 
as the angular distribution of galaxies on the sky (Coles & Plionis 1991; Gott et al. 1992), the 
distribution of galaxies in slices of the universe (Park et al. 1992; Colley 1997; Park, Gott, & Choi 
2001a), and the CMB temperature fluctuation field (Coles 1988; Gott et al. 1990; Smoot et al. 
1994; Kogut et al. 1996c; Colley, Gott, & Park 1996; Park et al. 2001b). 

We use the two-dimensional genus statistic introduced by Gott et al. (1990) as a quantitative 
measure of the topology of the CMB anisotropy. For the CMB the genus is the number of hot spots 
minus the number of cold spots. Equivalently, the genus at a temperature threshold level v is 

a{y) = ^J Kds, (4) 

where k is the signed curvature of the iso-temperature contours C, and the threshold level v is the 
number of standard deviations from the mean. At a given threshold level, we measure the genus 
by integrating the curvature along iso-temperature contours. The curvature is positive (negative) 
if the interior of a contour has higher (lower) temperature than the specified threshold level. 

The genus of a two-dimensional random-phase Gaussian field is g{y) oc ve~ u ' 2 (Gott et al. 
1990). Non-Gaussian features in the CMB anisotropy will appear as deviations of the genus curve 
from this relation (P98; Park et al. 2001a, 2001b). Non-Gaussianity can shift the observed genus 
curve to the left (toward negative thresholds) or right near the mean threshold level. The shift Av 
of the genus curve with respect to the Gaussian relation is measured by minimizing the x 2 between 
the data and the fitting function 

G = A s v'e- U ' 2 / 2 , (5) 

where v' = v — Av, and the fitting is performed over the range —1.0 < v < 1.0. Non-Gaussianity 
can also alter the amplitudes of the genus curve at positive and negative levels differently causing 
\g{y ~ —1)| ^ \g(y ~ +1)|- The asymmetry parameter is defined as 

Ag = A H - A c , (6) 

where 



A H = / 



rv2 
9obsdv/ / g&tdv, (7) 

JU1 
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and likewise for Ac- The integration is limited to —2.4 < v < —0.4 for Ac and to 0.4 < v < 2.4 
for Ah- The overall amplitude A of the best-fit Gaussian genus curve g^ is found from x 2 -fitting 
over the range —2.4 < v < 2.4. Positive Ag means that more hot spots are present than cold spots. 
This kind of topological information can not be drawn simply from the one-point distribution. In 
summary, for a given genus curve we measure the best-fit amplitude A, the shift parameter Ais, 
and the asymmetry parameter Ag. 



5.2. Effects of Foregrounds on CMB Anisotropy Topology 

We use our mock foreground-contaminated CMB anisotropy maps for the MAP experiment to 
study the sensitivity of the genus statistic to non-Gaussianity of cosmological origin. To compute 
the genus we stereographically project the maps on to a plane. This conformal mapping locally 
preserves shapes of structures. Additional Gaussian smoothings are applied during the projection 
so that the projected maps have total smoothing scales of 0?3, 0?5 and 1?0 FWHM. We then exclude 
regions with \b\ < 20° to reduce the effect of Galactic emission. We present the genus curves as a 
function of the area fraction threshold level va, which is defined to be the temperature threshold 
level at which the corresponding iso-temperature contour encloses a fraction of the survey area 
equal to that at the temperature threshold level ua for a Gaussian field (Gott et al. 1990). The 
ua = level corresponds to the median temperature because this threshold level divides the map 
into high and low regions of equal area. All our analyses are based on genus curves measured as a 
function of ua- 

Figure 7 shows genus curves measured from the mock MAP CMB maps in the ACDM model 
with various foreground contaminations at 0?3 FWHM resolution (open symbols). The genus curve 
of the foreground-free map (noise-added; filled circles) is shown for comparison. The effect of 
instrumental noise on the genus is significant on both 0?3 and 0?5 scales (filled circles in Fig. 7a; 
also see the genus amplitude parameter A values in Table 3 described below). It should also be 
noted that the 0?3 scale genus is reduced by the finite pixel size of the MAP simulation (Melott et 
al. 1989). 

Table 3 lists the genus-related statistics A, Au, and Ag measured from CMB maps with 0?3, 
0?5, and 1?0 FWHM resolutions using eqs. (5) to (7). The uncertainties are la errors derived 
from the variances over the eight octants and account for cosmic variance. Here Avq an d Ago are 
the values from the noise- and foreground-free pure CMB anisotropy maps. Hence Av — Av§ and 
Ag — Ag§ are the net effects due to instrumental noise or foreground emission. Table 3 shows 
that Galactic emission significantly affects the shift parameter on all angular scales, at a level 
exceeding 2 a. However, their effect become much less serious when only the Galactic emission 
model uncertainty AG needs to be accounted for. FIR extragalactic source emission has minimal 
effect on the genus, on all angular scales. Radio point source emission results in positive genus 
asymmetry, particularly on intermediate angular scales (0?5 FWHM) at a statistical significance 
level of about 2 a. We have also measured the genus from the radio source-contaminated CMB 
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anisotropy map with 0?5 FWHM resolution after broadly excising pixels above 4.7 times the rms 
level of the noise-added CMB anisotropy map, and find the asymmetry still exceeds cosmic variance 
(Ag- Ago = 0.035 ± 0.020). 

When Galactic diffuse emission and bright radio point source emission are removed from the 
observed temperature fluctuation map using accurate foreground models, one may conclude that 
non-Gaussianity of cosmological origin is present if the measured genus-related parameters |Az/| > 
0.02 (0.04) or \Ag\ > 0.03 (0.08) at 0?3 (1?0) FWHM scale at 2 a. 



6. Conclusions 

We study the effects of foreground contamination and instrumental noise on the simulated 
MAP data through the correlation function and genus statistics. The foregrounds we consider are 
diffuse Galactic dust and free-free emissions, far infrared extragalactic sources, and radio point 
sources. The mock observations are made in a spatially-flat ACDM model. 

Our study shows that FIR extragalactic sources in the IRAS 1.2 Jy galaxy catalog and the 
COBE-DIRBFi map have negligible effect on the correlation function and the genus. Radio point 
sources also have little effect on the correlation function when sources with amplitudes larger than 
5 times the rms of the noise-added ACDM CMB anisotropy map are removed. However, on 0?5 
scales radio sources cause positive genus curve asymmetry, a non-Gaussian feature, at a statistical 
significance level of about 2 a compared to cosmic variance (Table 3). This is the scale where 
the smoothing scale is small enough that strong radio sources are resolved, and where it is large 
enough that instrumental noise is not a dominant anisotropy source. This non-Gaussian feature 
in the genus curve is not significantly reduced by removing sources above 5 a. We conclude that 
good radio point source data at 90 GHz is essential for an unambiguous detection of non-Gaussian 
behavior of the CMB anisotropy on this scale. 

Diffuse Galactic emission affects the correlation function much less than cosmic variance. But 
it can produce negative shifts of the genus curve, another non-Gaussian feature, which are 2 or 
3 times larger than cosmic variance. These effects are significantly reduced if model Galactic 
emissions (e.g., K96 or FDS99) are subtracted from the observed M^4Pdata, unless the amplitude 
of Galactic foreground emission in these models are in great error. The typical shift of the genus 
curve due to the uncertainty in the Galactic emission models is estimated from the difference map 
between the K96 and FDS99 models. On scales below 1?0 it is of order of the shift due to cosmic 
variance. Therefore, an accurate Galactic emission model is still necessary to make reliable tests 
of Gaussianity of the CMB anisotropy. With the present understanding of Galactic emission and 
radio sources, one will be able to detect non-Gaussianity of cosmological origin at the 2 a level if 
the shift parameter of the genus curve \Av\ > 0.02 (0.04), or the asymmetry parameter \Ag\ > 0.03 
(0.08) on 0?3 (1?0) FWHM scales. 

There are other kinds of foreground contamination emission mechanisms, in addition to those 
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discussed above, such as the thermal Sunyaev-Zel'dovich effect, the Ostriker-Vishniac effect, grav- 
itational lensing, and so on (see Refregier et al. 2000 for a review). Some of these will need to 
be dealt with if we are to fully utilize the cosmological information in the anticipated MAP and 
Planck data. 

We acknowledge valuable discussions with D. Finkbeiner, K. Ganga, Juhan Kim, Sungho Lee, 
and P. Mukherjee, and helpful comments from the anonymous referee. CGP and CP acknowledge 
support from the BK21 program of the Korean Government and the Basic Research Program of 
the Korea Science &; Engineering Foundation (grant no. 1999-2-113-001-5). BR acknowledges 
support from NSF CAREER grant AST-9875031. We also acknowledge support from NSF grant 
AST-9900772 to J.R. Gott. 
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Table 1. Contributions to Mock 90 GHz MAP Data at 0?3 FWHM Resolution 



Maps RMS (//K) 

CMB 100.3 

MAP noise 35.8 

CMB+noise 106.5 

Galactic emission (K96) 14.0 

Galactic emission (FDS99) 12.5 
AG (K96-FDS99) 6.1 

Radio sources 26.8 

Radio sources (> 5 a rejected) 16.7 
FIR extragalactic sources 2.0 



Table 2. Correlation Function Acoustic Valley Position Measured from Mock ACDM Skies 



Map 



No cut CMB 1?122 

+noise 1?124 

J6J > 20° +noise 1?113 

+noise+Galactic emission (K96) 1?109 

+noise+Galactic emission (FDS99) 1?109 

+noise+AG 1?114 

+noise+FIR extragalactic sources 1?114 

+noise+radio sources 1?112 
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Table 3. Genus Statistics Measured from Mock MAP Data 



Res. 


Maps 




A 


Au - Au 


Ag - Ag 


0?3 


CMB 




961 ± 23 


0.000 ± 0.007 


0.000 ±0.011 




+noise 




1746 ± 36 


0.001 ± 0.009 


-0.008 ±0.013 




+noise+Galactic emission 


(K96) 


1730 ± 32 


-0.017 ±0.009 


-0.017 ±0.015 




+noise+Galactic emission 


(FDS99) 


1734 ± 33 


-0.024 ±0.009 


-0.017 ±0.014 




+noise+AG 




1739 ± 31 


0.006 ± 0.008 


-0.003 ±0.013 




+noise+radio sources 




1775 ± 123 


-0.001 ±0.010 


0.003 ± 0.013 




±noise±FIR extragalactic 


sources 


1747 ± 36 


0.002 ± 0.009 


-0.007 ±0.013 


0?5 


CMB 




553 ± 18 


0.000 ±0.012 


0.000 ± 0.020 




+noise 




614 ± 18 


0.001 ±0.012 


0.003 ± 0.021 




+noise+Galactic emission 


(K96) 


605 ± 15 


-0.025 ±0.013 


0.000 ± 0.024 




+noise+Galactic emission 


(FDS99) 


606 ± 17 


-0.034 ±0.013 


0.003 ± 0.022 




+noise+AG 




612 ± 16 


0.010 ±0.012 


0.002 ± 0.022 




+noise+radio sources 




629 ± 49 


-0.011 ±0.017 


0.041 ±0.020 




+noise+FIR extragalactic 


sources 


615 ±18 


0.000 ± 0.012 


0.004 ± 0.021 


1?0 


CMB 




181 ± 10 


0.000 ± 0.019 


0.000 ± 0.036 




+noise 




184 ± 10 


0.001 ±0.019 


0.003 ± 0.037 




+noise+Galactic emission 


(K96) 


180 ±8 


-0.033 ±0.020 


0.008 ± 0.041 




+noise+Galactic emission 


(FDS99) 


182 ±9 


-0.044 ±0.019 


0.001 ± 0.040 




+noise+AG 




183 ±9 


0.010 ±0.019 


0.025 ± 0.039 




+noise+radio sources 




189 ± 18 


-0.004 ± 0.029 


0.034 ± 0.045 




+noise+FIR extragalactic 


sources 


184 ± 10 


0.002 ± 0.019 


0.004 ± 0.037 
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Fig. 1. — FIR-to-microwave Galactic emission spectra at \b\ > 20°. K96 Galactic emission spec- 
trum (thick solid curve) is the sum of dust (dotted curve) and free-free (dashed curve) emission 
components, inferred by cross-correlating the DMR and DIRBE maps. Upper and lower thin solid 
curves are upper and lower limits of the K96 emission model. The thick long-dashed curve is the 
FDS99 dust spectrum, derived by averaging all the temperature-varying dust spectra at each line 
of sight and normalized with K96 spectrum at 3000 GHz (100 fim) for direct comparison. 



21 



0.5 



-0.5 - 




0.5 



-0.5 - 



Fig. 2. — Stereographic projections of 289 DIRBE FIR sources found in (a) north (152 groups) and 
(b) south (137 groups) hemispheres. The x and y coordinates are defined by x = cos I tan(7r/2 — 6)/2, 
y = sin / tan(7r/2 — b)/2 for the north hemisphere, x = cos 1/ tan(-7r/2 — b)/2, y = — sin // tan(-7r/2 — 
b)/2 for the south hemisphere. The three circles in each panel, from outside to the center, indicate 
|6| = 0°, 30°, 60°, respectively. The Galactic center is located at (x,y) = (1,0) in both panels. 
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Fig. 3. — The ACDM model foreground-contaminated CMB anisotropy temperature fluctuation 
mock map for the 90 GHz MAP channel, smoothed to 0?3 FWHM angular resolution. Foreground 
contamination here is modelled by the K96 Galactic emission model. 
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Fig. 4. — (a) Two-point angular correlation functions measured from the noise-free (filled circles, 
solid curve) and noise-added (open circles, dotted curve) ACDM temperature anisotropy maps at 
0?3 FWHM resolution. The theoretical CMB anisotropy angular correlation function for the ACDM 
model (dot-dashed curve) is shown for comparison, (b) Noise-added ACDM temperature two-point 
angular correlation functions showing the effect of different Galactic plane cuts. For |6| > 40° 
(filled triangles, dot-dashed curve), the angular correlation appears to be strongly biased by sample 
variance. 
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Fig. 5. — Two-point angular correlation functions measured from the noise-added ACDM model 
CMB anisotropy maps contaminated by (a) Galactic emission, (b) the Galactic emission model un- 
certainty, (c) FIR extragalactic source emission, and (d) radio point source emission. All correlation 
functions are measured from the Galactic plane cut maps (\b\ > 20°). In all panels, the foreground- 
free correlation function (solid curve) measured from the noise-added ACDM model map is shown 
for comparison. Note that 5 a rejection of bright radio point sources effectively restores the radio 
source- free result (d, open circles, dot-dashed curve). 
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Fig. 6. — Angular correlation functions in the region of acoustic valleys. Error bars are la uncer- 
tainties estimated from 15 independent correlation function measurements with different random 
seeds to pick up a million pairs of points differently. Solid, dotted, and dashed curves are 4 th -order 
polynomial fits used to measure acoustic valley positions, listed as 9 V in the lower-left part of each 
panel. Error bars for FDS99 case in (b) are omitted for simplicity. 
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Fig. 7. — Two-dimensional genus per steradian measured from the noise-added ACDM model 
CMB anisotropy map (0?3 FWHM, \b\ > 20°) and from map(s) (a) without noise and foreground 
contamination (open circles), (b) with noise and Galactic foreground emission (open circles for 
the K96 model, open triangles for the FDS99 model), (c) with noise and the Galactic emission 
model uncertainty (open circles), and (d) with noise and radio point source emission (open circles). 
Filled circles in each panel are the genus measured from the noise-added CMB map, with error 
bars derived from the variance of the genus computed independently in each of the eight octants 
of the sphere. Dotted curve in each panel shows the functional form expected for a random-phase 
Gaussian field, Nve~ u ' 2 , and has been fitted to the filled circles by adjusting N. Small panels in 
the upper-left region of the large panels from (b) to (d) show genus deviations of the foreground- 
added CMB maps relative to the noise-added CMB result. Genus of the CMB map contaminated 
by FIR extragalactic source emission is almost the same as the foreground-free result and is not 
shown here (see Table 3). 



